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A central ingredient in establishing the properties of the newly discovered Higgs-like boson is to 
isolate its production via vector boson fusion (VBF). With the typical experimental selection cuts, 
the VBF sample is contaminated by a ~ 25% fraction from Higgs + 2 jet production via gluon fusion 
(ggF) which has large perturbative uncertainties. We perform a detailed study of the perturbative 
uncertainties in the NLO predictions for pp H + 2 jets via ggF used by the ATLAS and CMS 
collaborations, with the VBF selection cuts of their current _ff — >■ 77 analyses. We discuss in detail 
the application of the so-called "ST method" for estimating fixed-order perturbative uncertainties 
in this case, and also consider generalizations of it. Qualitatively, our results apply equally to other 
decay channels with similar VBF selection cuts. Typical VBF selections include indirect restrictions 
or explicit vetoes on additional jet activity, primarily to reduce non-Higgs backgrounds. We find 
that such restrictions have to be chosen carefully and are not necessarily beneficial for the purpose of 
distinguishing between the VBF and ggF production modes, since a modest reduction in the relative 
ggF contamination can be easily overwhelmed by its quickly increasing perturbative uncertainties. 



I. INTRODUCTION 



With the discovery of a Higgs-hke boson by the ATLAS 
and CMS collaborations [l|, 0]; a central ingredient in 
measuring the properties of the new particle is to separate 
out the different production mechanisms via gluon-gluon 
fusion (ggF) and vector boson fusion (VBF). 

To maximize the signal sensitivity, the experimental 
analyses separate the data into various exclusive selection 
categories, based on the number of jets ("jet bins") and 
other criteria. One selection category, designed to isolate 
a clean VBF signal, is the production in connection with 
two jets that are widely separated in rapidity. To fur- 
ther enhance the VBF signal and reduce non-Higgs back- 
grounds as well as the sizable contamination from ggF 
production, additional kinematic selection criteria are ap- 
plied. A characteristic feature of the VBF process is that 
it is accompanied by few extra gluon emissions, because 
of its color structure and incoming quarks. The same 
is not the case for ggF production or generic non-Higgs 
backgrounds, so the VBF signal tends to be most signifi- 
cant in the exclusive 2-jet region of phase space with two 
forward jets and little additional radiation. Therefore, 
when optimizing its significance, whether in a cut-based 
approach or via multivariate techniques, one dominantly 
selects events from this region. 

In general, placing a restriction on additional real emis- 
sions induces Sudakov logarithms at each order in the 
perturbative series. In the limit of very tight restrictions 
the logarithms become large and must be resummed to 
all orders in the strong coupling constant to obtain 
a meaningful perturbative prediction. Often, the experi- 
mentally relevant region is an intermediate one where the 
logarithmic corrections are already sizable but their re- 
summation is not yet strictly necessary and a fixed-order 
expansion can still be applied. In this case, however, it is 
important that the possible sizable effects of higher-order 
logarithms are reflected in the perturbative uncertainty 



estimate for the fixed-order prediction. 

Due to the incoming gluons in ggF and the associated 
large color factor, the logarithmic corrections in this in- 
termediate region are indeed sizable. This was shown 
explicitly for the 0-jet bin in Refs. 0, 0] which has been 
computed and studied extensively to NNLO [5l-[To|. A 
numerically important ggF contribution to the VBF-like 
2-jet selection is the partonic gg — >■ Hgg process, where 
both incoming and outgoing gluons generate additional 
radiation. When restricting that radiation by forcing the 
kinematics into the exclusive 2-jet region, the logarithmic 
corrections can be expected to be at least as large and 
likely larger than in the 0-jet case. Hence, the perturba- 
tive uncertainties have to estimated carefully, in partic- 
ular since here NNLO predictions are not available. 

In Ref. Q a simple method was devised that explicitly 
takes into account the size of the logarithmic corrections 
in the fixed-order uncertainty estimate (for which a sim- 
ple scale variation in the exclusive jet cross section is 
insufficient). It has been adopted in Refs. [HI [13, and is 
being employed in various exclusive analyses at the LHC 
and the Tevatron. It is sometimes referred to as "ST 
method" . Here we apply this method to provide robust 
uncertainty estimates for the NLO calculation [l^ [13] 
in the exclusive 2-jet bin in ggF production that is cur- 
rently being used in the ATLAS and CMS Higgs analyses. 
An independent NLO calculation has been performed re- 
cently in Ref. ISi]. We also discuss the application to 
more general cuts restricting to the exclusive 2-jet region. 
The application to multivariate selection techniques will 
be discussed in a forthcoming publication [l6| . 

In principle, by performing a higher-order logarith- 
mic resummation one can gain additional information, 
which allows for refined perturbativepredictions and un- 
certainty estimates, see e.g. Refs. [3, [T7l - l22j for recent 
applications to Higgs plus or 1 jets. Although techni- 
cally more demanding, similar methods could be used in 
the future to provide improved predictions for the exclu- 
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sive gg ^ H + 2 jets cross section. An important step 
in this direction was made in Refs. [1^, [l^, where the 
fixed NLO predictions {oi gg ^ H + 2 jets are matched 
to a parton shower, which on top of the NLO corrections 
provides a leading-logarithmic (LL) resummation in the 
exclusive 2-jet region. Here, our NLO uncertainty analy- 
sis provides an important baseline for future studies, as it 
is often difficult to obtain a robust resummed uncertainty 
estimate from a LL resummation alone. 

In the next section, we give a general discussion of 
jet binning uncertainties, reviewing and extending the 
method to estimate fixed-order uncertainties introduced 
in Ref. [Ij. In Sec. IIIII we discuss and validate its ap- 
plication to the ggF contribution in the VBF-like 2-jet 
selection. In Sec. IIVI we present our results for the NLO 
uncertainty estimates, implementing the VBF selection 
cuts used in the iJ ^ 77 analyses by the ATLAS [H HI] 
and CMS [27[| experiments. Our findings apply equally 
to other decay channels where similar VBF selection cuts 
are applied. In Sec. HV] we also discuss the implications 
of the ggF uncertainties for the VBF-ggF separation. We 
conclude in Sec. |Vl 



II. JET BINNING UNCERTAINTIES 



with pn+1 = for a tree-level 7V-parton state. Typi- 
cal examples would be the pt of the N + 1st jet or the 
total \pt\ of the underlying A^-jet system. The two jet 
bins then correspond to the integrals of the differential 
spectrum da/dpN+i above and below some cut, 



cr>N = 



dp 



N+1 -. \- / dpN+i - 

clpAT+i Jpcut apN+1 

^7w(p™*)+'T>Ar+i(p-"). (2) 



In general, the jet bin boundary could be a much more 
complicated function of phase space, for example in a 
multivariate analysis. In Sec. IIVI we will also consider the 
next-to-simplest case of a two-dimensional rectangular 
cut on two kinematic variables. 

We are interested in the uncertainties involved in the 
binning. The covariance matrix for {ctat, (T>Ar+i} is a 
symmetric 2x2 matrix with three independent param- 
eters. A convenient and general parametrization is to 
write it in terms of two components, 



C = 



(A?v) 



^Af ^>Af-|-l 



A2 , 

cut 



cut 



(3) 



Consider the inclusive A^-jet cross section, a>N, for 
some process containing at least N jets. We will assume 
that cr> AT is a sufficiently inclusive quantity such that it 
can be computed in fixed-order perturbation theory. We 
are interested in the case where a>N is divided up into 
a corresponding exclusive A^-jet cross section, cr^v, and a 
remainder a>N+i, 

(T>N — cr Ar(excl. cut) -I- cr> AT+i (inverse excl. cut) . (1) 

All three cross sections here have the same selection cuts 
applied that identify the leading A^ signal jets. What de- 
fines ctat to be "exclusive" is that the additional exclusive 
cut applied to it restricts the phase space of additional 
emissions in such a way that is dominated by con- 
figurations close to the Af-parton Born kinematics. In 
particular, at leading order (LO) in perturbation theory 
^>Ar = '^N^^ while relative to these CT>Ar+i is suppressed 
by 0{as)- In other words, ct>jv+i requires at least one 
additional emission to be nonvanishing. Hence, we can 
consider it an inclusive (A^ + l)-jct cross section with at 
least A^ -I- 1 jets.^ 

In the simplest case, iT>Ar is divided into the two jet 
bins ctat and a>isr+i by using a cut on some kinematic 
variable, pn+i, which characterizes additional emissions. 



^ It should be stressed that o'>jv+i here is defined by inverting the 
exclusive cut that defines apf and so does not necessarily require 
the explicit identification of another well-separated jet via a jet 
algorithm. The variables PTHjj and A{j>H-jj we will use are 
examples of this. 



Here, the first term is an absolute "yield" uncertainty, 
denoted with a superscript "y" , which (by definition) is 
100% correlated between the two bins and o'>Ar+i. 
The second term is a "migration" uncertainty between 
the bins and corresponds to the uncertainty introduced 
by the binning cut. It has the same absolute size, Acut, 
for both bins and is 100% anticorrelated between them, 
such that it drops out when the two bins are added. 
Hence, the total uncertainty for each bin is given by 



(A]; ^2 



(4) 



while the total uncertainty on their sum, i.e. on (t>n, is 
given by the total yield uncertainty. 



A^ 



(5) 



Considering the perturbative uncertainties, the basic 
question is how each of the uncertainties in Eq. ^ 
can be evaluated. The fixed-order prediction provides 
us with two independent pieces of information, namely 
the variations obtained by the standard scale variations, 
which we denote as A>^, A^, A^j^+ii ^tnd which sat- 
isfy A>^ = A^ -I- A>^_,_^. Following Ref. we start 
by assuming that the standard fixed-order scale varia- 
tions can be used to obtain a reliable estimate of the 
total uncertainties in the inclusive cross sections (which 
is of course the common assumption underlying any in- 
clusive fixed-order calculation). Hence, we impose the 
two well-motivated boundary conditions 



A>Ar = A'^ 



>N ' 



(6) 
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Together with Eqs. (jH) and ([5]) these lead to 



i) 
ii) 



A^ 



(A^ 



>N+1) 



(7) 



Thus, the question is how to divide up A>^^-^ between 
^>Af+i ^^'^ ^cut in order to satisfy condition ii). Con- 
dition i) then determines A^. The nontrivial effect Acut 
can have is on the size of An as weh as on the off-diagonal 
entries in Eq. which determine the correlation be- 
tween A^v and A>jv+i. 

Clearly, the simplest is to neglect the effect of Acut 
altogether and to directly use common scale variations 
to estimate the uncertainties, i.e., to take 



A^ 



AM 



AM 



AM 



A^ 



AM 



Acut = , 
which leads to 

direct: C 



/AM \2 AM AM 



N 

A^ir Af" 



N '-^>N+1 



(8) 



(9) 



Note that since a>N+i starts at higher order in per- 
turbation theory than cr> n , its relative uncertainty 
'-^>N+i/'^>N+i will typically be (much) larger than 
CT>jv's relative uncertainty Ay^/a>M- This means one 
cannot simply apply the latter as the relative yield un- 
certainty in each bin by taking A\ = (A>^/cr>Ar)i7i, as 
this would violate the condition A>jv+i = A^^_|_-^. This 
point has already been emphasized in earlier studies Q. 

The direct scale variation choice is reasonable as long 
as the effect of Acut is indeed negligible. It is certainly 



where as denotes the Born cross section and L is a Su- 
dakov logarithm, e.g. for Eq. ^ L = ln(j)'^"'/(5) where 
Q ^ rriH is a typical hard scale. As the binning cut 
becomes tighter {p'^^^ becomes smaller) the logarithms 
grow in size. Once the logarithms are 0(1) numbers one 
is in the transition region and the logarithms will start 
to dominate the perturbative series of a>N+i, and there 
will be sizable cancellations in a^f between the perturba- 
tive series for a>N and the logarithmic series in aypf+i. 
Eventually the logarithms will grow large enough to over- 
come the as suppression and (Tn becomes negative, at 
which point one is in the resummation region and the 
fixed-order expansion has broken down. 

The perturbative migration uncertainty Acut can be 
directly associated with the perturbative uncertainty in 
the logarithmic series induced by the binning, and so 
should not be neglected once the logarithms have a no- 
ticeable effect. In particular, as demonstrated in Ref. Q, 
the simple choice in Eqs. (jH) and ([9|) can easily lead to 
an underestimate of Ajv in the region where there are 
sizable numerical cancellations between the two series in 
cr>Ar and (t>jv+i- Since in this region the dominant con- 
tribution to (T>N+i comes from the logarithmic series, 
varying the scales in (t>7v+i directly tracks the size of 
the logarithms, which means we can use Acut = A^j^^^ 
as an estimate for the binning uncertainty, which is the 
basic idea of Ref. 3 . From Eq. ([7]) , we then find 



A^ - A^ 

A — A'^ 
^cut — ^>7v+l 



A'' - n 



(11) 



such that 



justified if numerically A>^ ^ A^^_|_-^, since any uncer- 
tainty due to migration effects can be at most as large as 
^>N+i i^y virtue of condition ii). This can happen for ST: C 
example when A^j^ is sizable due to large perturbative 
corrections in a>N and/or the binning cut is very loose 
(i.e., is cutting out only a small fraction of phase space) 
such that a>N+i is numerically small to begin with. 

In perturbation theory, the effect of the binning cut is 
to introduce Sudakov double logarithms in the perturba- 
tive series of and (t>7v+i, which have opposite sign 
and cancel in the sum of the two bins, schematically 



(A^ 



(12) 



^ is now used as Acut, the effective outcome 



Since A^^j^^ 

is that one treats A^^ and A^j^j_^_-^ as uncorrelated. 

More generally, we can introduce a parameter < 
p < 1, which controls the fraction of A'^j^^^ assigned 
to A^^^j^, such that 



<7>N ~aB[l + as+al+ ©(a^)] , 
<y>N+i — (Tb [as{L^ + L+l) 

+ al{L'^ + L^ + L^ + L + l) + ©(a^L^)] , 
(Jn = a>N — cr>N+i , (10) 



A^ - A*" - n A^ 



Acut = Vl A'^jy_^^^ , 



which leads to 



A^ -nAf^ 



>N+1 ' 



(13) 



ST (p): 



C 



(A^ 



(PA'>^-A^t 



AM \AM 



>7V+1 



^A^ 



(14) 



4 





ATLAS 


CMS loose 


CMS tight 




anti-fcr R = 0.4 


anti-ZcT R = 0.5 


anti-fcr R = 0.5 


2-jet selection 


PTj >25GeV for |r7j |<2.5 


jet 1: pTj> 30 GeV, \Vj\<'i-7 


PTj >30GeV, |77i|<4.7 




PTj > 30 GeV for 2.5 < | < 4.5 


jet 2: pTj>20GeV, |7?i|<4.7 




^Vjj = I%1 - ^J2l 


> 2.8 


> 3.0 


> 3.0 




> 400 GeV 


> 250 GeV 


> 500 GeV 


\VH - ivn +'7j2)/2| 




< 2.5 


< 2.5 


A(j}H-jj 


> 2.6 


> 2.6 


> 2.6 



TABLE L VBF selection cuts we use, corresponding to the H ^ -j-j analyses by ATLAS and CMS CMS loose 

excludes events that pass CMS tight. The cut on A(f)H-jj in the last row is treated special as an exclusive binning cut. 



From this one can easily see that p corresponds the cor- 
relation between A^jy and A^jy^j^. The choice p = 1 
would be equivalent to the case in Eq. ([9]), while p — Q 
reproduces Eqs. pT|) and p^ . Hence, from the above 
arguments one should take p to be small. In the next 
section, we will explore the dependence on p in the ST 
method. We will see that all choices p ^ 0.4 give very 
similar results, so for our results in Sec. IIVI we will use 
the default choice p — Q. 

As a final comment, note that in general one could also 
take /9 to be a function of the binning cut. For example, 
at large p™'^ the logarithms become small, in which case 
one might want to reproduce the direct scale variation 
uncertainties in Eq. However, in this limit typically 
A>jY+i becomes much smaller than A^jy: which makes 
the precise choice of p irrelevant there, and so it is con- 
sistent to use a fixed p = everywhere. 



III. APPLICATION TO gg ^ H + 2 JETS 

We now discuss the application of our method to the 
case oipp — > H + 2 jet production via gluon fusion (which 
for simplicity we denote as gg ^ H + 2j, where a sum 
over all possible partonic channels is implied). We will 
study the uncertainties in the exclusive -I- 2 jet cross 
section as a function of two kinematic variables, pxHjj 
and A(j>H-jj- 

We take Vs = 8TeV and mn = 125 GeV. We use 
MCFM [H [H, il to compute the NLO cross sec- 
tion, with the ggH effective vertex in the infinite top 
mass limit. We then rescale the cross section with the 
exact rrit dependence of the total Born cross section, 
crB(mt)/crs(oo) = 1.0668. We use the MSTW2008 [29^ 
NLO PDFs with their corresponding value of as{mz) = 
0.12018. For all our central value predictions we use 
jjLr = Pf = tuh, which was also used in Refs. [isl . [3]. 
The scale variations in the inclusive cross sections are 
discussed below in Sec. IHIBI 

In our analysis we implement the 2-jet selection and 
VBF selection cuts summarized in Table HI which are 
taken from the current ATLAS and CMS 77 anal- 

yses. However, note that we consider the cross section 



for the production of an on-shcU Higgs boson, without 
including any branching ratios or cuts on the Higgs de- 
cay products. 

A. Variables 

1- PTHjj 

We define PTHjj as the magnitude of the total trans- 
verse momentum of the Higgs-dijet system, 

PTHjj = \PT jl + PTj2 + PTh \ ■ (15) 

At Born level, PTHjj — and so applying a cut pthjj < 
p'thj] restricts the phase space to the exclusive 2-jet re- 
gion. At NLO PTHjj is equivalent to the px of the third 
jet, so it is a useful reference variable for a pr-veto on ad- 
ditional emissions, such as the central jet vetoes applied 
in the H — >■ WW and H ^ tt VBF analyses (see e.g. 
Refs. (S^-Ull)-^ It is also considered directly for example 
in the latest H — >■ rr analysis [3l[ . 

The exclusive 2-jet cross section (T2 [pTHjj < p'thjj ) i^ 
shown in the left panel of Fig. [1] as a function of p'rikjj 
and using three different combinations of the factoriza- 
tion and renormalization scales pr and pf. The solid 
line and blue band correspond to pr = itih and varying 
Pf — {2, 1, l/2}m/f. Similarly, we vary pf while keep- 
ing Pr = mn/^ for the dark green band and pr = 2m// 
for the light green band. One can see that the biggest 
variation is due to the pr variation, while the pj vari- 
ation only has a subdominant effect, which was already 
noticed in Ref. [l3| • Therefore, for simplicity we will take 
Pr ~ Pf = p and vary p = {2, 1, l/2}mH when showing 
the direct scale variations as reference in the following. 



^ The central jet veto is applied to reconstructed jets at central 
rapidities, which at low px values can be heavily influenced by 
underlying event, pile-up, and detector effects. Since none of 
these effects can be accounted for by the NLO calculation, we did 
not attempt to study an explicit central jet veto here. Instead, 
we concentrate on PTHjj, which is cleaner as it only requires 
information about the two signal jets and the Higgs candidate. 
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FIG. 1: Exclusive 2-jet cross section using the ATLAS VBF selection for various scale choices as a function of PxHjj (left panel) 
and TT — A(j>'H-jj (right panel). 



We write the exclusive 2-jet bin defined by this eut in 
terms of the inclusive 2-jet cross section, (t>2, and the 
inclusive 3-jet cross section with the cut inverted as, 

<^2{pTHjj < Prhjj) = cr>2 - <y>3{PTHjj > PTHn) ' (1^) 

where in all cases the remaining VBF selection cuts in 
Table U are applied (excluding the cut on A(j)H^j in this 
case). 

The restriction on pxHjj is infrared sensitive and in- 
duces Sudakov logarithms of the form L = ln{p^^j^ / mu) 
in the perturbative series of (T2 and cr>3. In Fig. [1] we 
see that the veto starts to have a noticeable effect be- 
low PTHjj ^ 50 GeV, where the different scale variations 
start crossing and we start to see cancellations between 
cr>2 and a>3. In the region below PTHjj ^ 20 GeV the 
logarithms have grown large enough for the NLO cross 
section to go negative and the fixed-order perturbative 
expansion to break down. In the intermediate region in 
between, the fixed-order prediction can still be used, but 
the direct scale variation does not provide a reliable un- 
certainty estimate as it does not properly take into ac- 
count the effect of the binning cut. 



so the constraint A(pH-jj > ^(pW-jj forces the kinemat- 
ics into the exclusive 2-jet region and restricts additional 
emissions. Hence, it behaves similar to p'rHjj ^^"^ f^'" 
TT — A(f>'j!^^jj — induces large logarithms in the pertur- 
bative series. The exclusive 2-jet cross section in terms 
of A(ff^'^^j is written as 

a2{l^(t)H-n > ^<t>H-n) = ^>2 - a>3{A(j,H-n < , 

(18) 

with the remaining VBF cuts applied in all three cross 
sections. The right panel of Fig. [1] shows (T2{A(j)H-jj > 
A05i-j j ) plotted as a function of tt — A(/)^_*j-j , where one 
can clearly see the very similar behavior to the p'rHjj 
case in the left panel. Here, the exclusive cut on A(f>H-jj 
starts having a noticeable effect below tt — Aipn-jj ^ 0.6 
and the fixed-order perturbative expansion breaks down 
below around tt — ^ 0.2. In the transition region in 
between the direct scale variations again do not provide 
a meaningful uncertainty estimate, because they neglect 
the effect of Acut- 



B. Inclusive Scale Uncertainties 



2. A(t)H-jj The two fixed-order scale variation uncertainties we 

require as inputs are A>2 and A^^.3. In Fig. [H one can 
As shown in Table H the VBF category in the if 77 already see that the scale variationls asymmetric at large 

values of and tt — Acff^z.^^ In Fig. we show 

the scale dependence of the inclusive 2-jet cross section, 
CT>2, where we plot it over a range of 1/4 < iir/mn < 4 
for three different values of fif. We take /i/ = /i,. = 
TO/f , corresponding to the fXr/mn — 1 point on the blue 
solid line, as our central value for (J>2, and consider the 
range 0.5 < /Zr/m_y < 2 to estimate the inclusive scale 
uncertainty. The maximum deviation from the central 
value is given by the green dotted curve for /iy = /i^ = 
2mH- We use this maximum variation to construct a 



analyses by ATLAS and CMS includes a cut A4)H-jj > 
2.6 radians (150 deg) (where the Higgs momentum is rep- 
resented by the total momentum of the diphoton system) . 
Taking the beam direction along the z-axes, Acjjn-jj is 
defined as 



cos AcpH-^jj = — 



(PTjl +PTj2) -PTH 



\PTjl +PTj2\\PTH\ 



(17) 



Momentum conservation in the transverse plane implies 
that events with only two jets always have A(t>H-jj ~ tt, symmetric uncertainty A>2: as shown by the uncertainty 
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FIG. 2: Inclusive 2-jet cross section over a range of ^r/mu for ATLAS VBF selection (left panel) and CMS loose selection 
(right panel). The three curves show different values of /i/. The blue solid, green dotted, and green dashed curves correspond to 
/i/ = niH, fJ-f = ^iriH, and = m_r//2, respectively. The uncertainty bars show the inclusive 2-jet scale variation uncertainty. 




FIG. 3: Inclusive 3-jet cross section as a function of p'thh (l^ft panel) and tt — Atjf'n^j (right panel) for the ATLAS VBF 
selection. The outer solid green lines show the inclusive 3-jet scale variation uncertainty after symmetrization. 



bar in the figure. It corresponds to a relative uncertainty 
at NLO of 21%, which is similar to what was found in 
earlier studies [13 where a somewhat looser VBF 
selection was used. The corresponding uncertainty at 
LO is +76% and -40%. 

In Fig. 131 we illustrate the scale variation uncertainties 
for the inclusive 3-jet cross section, (t>3, for both p'rHjj 
and Acjfj^^jj , and using the ATLAS selection as example. 
(The results for cr>3 with the CMS selections look very 
similar except for the different overall scale.) The blue 
line shows the cross section for fir ~ fJ-f — fnn , which we 
take as the central value for (T>3. The green dashed and 
dotted lines show the scale variations fir = fif = mn/'i 
and fir = fJif — 2m respectively. For simplicity, we 
symmetrize the uncertainty by taking half of the differ- 
ence between the up and down variations as the inclusive 
3-jet scale uncertainty A^^.^, i.e., we keep the size of the 
band and move it to be symmetric about the central blue 



line, which is shown by the outer solid green lines. The 
relative uncertainty is of 0(70%) and almost indepen- 
dent of p'rHjj ^'^'i ^^'H-jj ■ This rather large uncertainty 
is not too surprising, since this is a leading-order H + 3j 
cross section, which starts at a^. 

C. Exclusive Uncertainty 

Having obtained the perturbative uncertainties 
and in the inclusive cross sections from the usual 
scale variation, we now study the resulting uncertainty 
A2 in the exclusive 2-jet cross section according to the 
discussion in Sec. [ill From Eq. ()14|) we have in general 

Al ^ (A^^)' + (A^a)' - 2p A^,A^3 , (19) 

where p is the assumed correlation between A^^.2 and 
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FIG. 4: Perturbative uncertainties in the exclusive 2-jet cross section with the ATLAS VBF selection as a function of p"^^ 
(left panel) and vr — (right panel) for different choices of the correlation parameter p. Our default choice is p = Q. 



1. Dependence on p 

We first investigate the dependance on the choice of 
p. In Fig. 13] we show the uncertainty in the exclusive 2- 
jet cross section as a function of p'thh ^"^^ ^'/^'n-jj 
different values of p from to 0.95. The outermost solid 
curves show the uncertainty obtained with our default 
choice p = 0, which effectively assumes that A>2 and 
are uncorrclatcd. For p ^ 0.4 the results are not 
very sensitive to the precise value of p, which is reassuring 
and shows that p = is in fact a safe choice on the 
conservative side. 

As p increases further, the uncertainty bands in the 
transition region keep shrinking, and for p = 0.95, shown 
by the innermost dot-dashed lines, pinch near PxHjj — 
30 GeV and tt- A0^^^^ ~ 0.3. (For p = 1 the uncertainty 
goes exactly to zero around these points.) This is because 
for p — >■ 1, A>3 and become 100% correlated, which 
is equivalent to the case of direct scale variation. (The 
only difference compared to the direct scale variations 
we saw in Fig. [T] is that here we symmetrized the scale 
variations.) 

One can also see that for large cut values, where the 
veto is not relevant and we approach the inclusive 2-jet 
cross section, the choice of p becomes irrelevant, because 
the absolute size of A^g becomes numerically negligible 
compared to A^^.2. 



2. Comparison to Efficiency Method 

Another prescription to obtain fixed-order uncertainty 
estimates for exclusive jet cross section, which is based on 
using veto efficiencies, was applied in Ref. [l3| to the 0-jet 
case at NNLO. We will refer to it as "efficiency method" . 

In Ref. [S^l it was shown that for the case of H + jets 
at NNLO the ST method and efficiency method yield 



very similar uncertainties, providing a good cross check 
on both methods. 

The starting point in the efficiency method is to write 
the exclusive jet cross section in terms of the correspond- 
ing inclusive jet cross section times the corresponding 
exclusive efficiency, i.e., applied to our 2-jet case, 

(72 = 0->2 1 ^ = 0->2 X €2 , 

V cr>2/ 

CT>3 = <7>2 (1 - £2) , (20) 

where the logarithmic series induced by the jet bin ning 
now only affects the efficiency. The basic assumption 
one then makes is to treat the perturbative uncertainties 
in (T>2 and 62 as uncorrelated (which one can think of 
as a multiplicative version of the ST approach). One 
should be aware that this method does not satisfy one 
of our starting conditions, namely the total uncertainty 
A>3 for cr>3 will not be given by its standard scale vari- 
ation A^i.3 anymore. Nevertheless, it is a useful way to 
gain additional insights into the size of higher-order cor- 
rections. 

The 2-jet efficiency £2 = 1 — c>3/o'>2 is still an ex- 
clusive quantity. Similar cancellations between the two 
perturbative series for a>2 and cr>3 can happen in their 
ratio than in their difference, so the direct scale variation 
for £2 might not provide a reliable uncertainty estimate. 
To circumvent this, in Ref. [iTj the perturbative uncer- 
tainty in e is instead estimated by using three different 
schemes for how to write the perturbative result for e, 
which are all equivalent up to the desired order in as, 
but differ in the higher-order terms that are retained or 
not. 

The inclusive 2-jet and 3-jet cross sections have the 
following perturbative structure 

0'>2 = [o'>2 + (^>l + 0?s <^>2 + ^(Q^s)] ' 

a>3 = al + ai^ + 0{al)] . (21) 
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FIG. 5: Exclusive 2-jet efficiency for different schemes in the efficiency method for Pthh (left panel) and tt — Acjf^ljj (right 
panel) using the ATLAS VBF selection. 
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FIG. 6: Comparison of the ST method with the efficiency method for Prikjj (left panel) and tt — A^^^^j (right panel) using the 
ATLAS VBF selection. The exclusive scale uncertainties from both methods are consistent with each other. The uncertainties 
from the efficiency method are very close to those from the ST method with p = 0.4. 



At NLO, the pieces we have available are a^l, o->2' ^^'^ 

(T>3 . In scheme (a) one defines the efhciency by keeping 
the full expressions in numerator and denominator, which 
at NLO gives 



= 1 - 



0'>3 
CT>2 



AO) 



j. + Oiai). (22) 



.(1) 



In scheme (b) one keeps the same number of terms in the 
perturbative series in the denominator as in the numera- 
tor, which in our case amounts to dropping the a^l term 
in the denominator. 



(0) 

>3 



(23) 



Finally, in scheme (c) one strictly reexpands the ratio 
to a given order in as, which to 0{as) unfortunately 



yields the same result as scheme (b) . To produce another 
expression with differing higher order terms, the closest 
scheme (c) analog we can do is to keep the 0{a'^) cross 
term that comes from expanding the denominator, so 



Ac) 



^>2 



^>2 
/O) 
^>2 



(24) 



In Fig.[5]we show the result for £2 in the three schemes 
for both pf}j^^ and Acjy'j^'^^^ using the ATLAS VBF se- 
lection. The central lines show the results for fir ^ fJ-f ~ 
fi = mn, while the bands are obtained from varying 
/i = {2,l/2}m// in each scheme. At NLO the central 
values from the three schemes are quite close and still lie 
within the direct scale variation of scheme (a), so their 
difference does not provide a useful uncertainty estimate 
here. The direct scale variation in scheme (b) is very 
small and in scheme (c) abnormally large (which is very 
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similar to what was seen in Ref. [TtI). Hence, in the end 
the most reasonable choice to get an uncertainty estimate 
for 62 is to just use the direct scale variation in scheme 
(a). 

In Fig. |6] we compare the results of the ST and effi- 
ciency methods for the exclusive 2-jet cross section 
for both vfHjj and A05f ^-j- using the ATLAS VBF se- 
lection. The blue solid curve shows our usual NLO cen- 
tral value, which is equivalent to the central value from 
scheme (a). The light orange solid curves are the uncer- 
tainties obtained in the efficiency method by combining 
the scale uncertainties A>2 with the direct scale vari- 
ations in 62"^ treating both as uncorrelated. The dark 
orange solid curves show the default ST uncertainties for 
/9 = 0, which are somewhat larger. The dashed lines 
show the ST uncertainties for p — 0.4, which agree al- 
most perfectly with the efficiency method. This result 
is not surprising. Basically, to obtain the £2 scale un- 
certainty we vary the scales correlated in tT>2 and cr>3, 
which has the effect of reintroducing a certain amount of 
correlation between and A^^g when considering (T2, 
which is also what a nonzero value of p does. Overall, the 
good consistency between the various methods gives us 
confidence in the reliability of our uncertainty estimates. 



IV. RESULTS 

In this section, we present our results for the exclusive 
pp — >■ 7J -|- 2 jet cross section via ggF at NLO, taking ST 
with p = as our method of choice to estimate the per- 
turbative uncertainties. All our inputs are summarized 
at the beginning of SecHTH The ATLAS, CMS loose, and 
CMS tight VBF selection cuts we apply are summarized 
in Table m 



A. — >■ + 2 Jets Cross Section 

1- PTHjj and l^4>H-ji 

In Fig. [7] we plot the result for the exclusive 2-jet cross 
section as a function of p'rHjj and A^^J';^^ for the AT- 
LAS, CMS loose, and CMS tight VBF selections. In aU 
our cross section plots the blue central line shows the 
central- value prediction obtained from p, = uih-, while 
the orange solid lines show our uncertainty estimate. For 
reference, the green dashed and dotted curves show the 
direct scale variation for p — mu/'i, and p = 2mH, re- 
spectively. 

The overall picture is very similar for all three VBF 
selections and both binning variables. For large values 
of p'rHjj or TT — A(f>'j!^*^jj, the cross section a>3 that is 
cut away becomes small and so the effect of A'^"' is neg- 
ligible. In this limit the uncertainties reproduce those 
in the inclusive 2-jet cross section, which here are deter- 
mined by the p = 2mH variation (cf. Sec. IIII Bl) . On 



the other hand, in the transition region, once the ex- 
clusive cut starts to impact the cross section, the direct 
scale variations cannot be used any longer to estimate 
uncertainties, which is exhibited by the crossing of the 
lines. As explained in detail in the previous two sections, 
the reason is that the direct scale variation only gives an 
estimate of the yield uncertainties, which effectively as- 
sumes the scale variations in the inclusive cross sections 
to be 100% correlated (corresponding to p = 1). At the 
same time it neglects the migration uncertainty in the 
binning, which becomes important as the exclusive cut 
gets tighter. In the ST procedure, this effect is taken 
into account explicitly, which thus gives more robust un- 
certainties for all values of p'rHjj or ^'t^'n^jj ■ 

In Table [ll] we quote results for the cross sections 
and their percentage uncertainties for specific cuts. For 
A<j)H-jj we use the current experimental value A(j)H-jj > 
2.6. Compared to the 21% in the inclusive 2-jet cross sec- 
tion with VBF cuts {a>2), we see a moderate increase in 
the uncertainty in a2iA(j)H-jj > 2.6) to 26% for ATLAS 
and CMS tight, and 24% for CMS loose. For prHjj we 
use a representative value of Prujj < 30 GeV, for which 
the uncertainties increase substantially to 44% and 49% 
for ATLAS and CMS tight, and moderately to 28% for 
CMS loose. Note that for a fixed cut the uncertainties 
increase with a tighter VBF selection. This is also clearly 
visible in Fig. [7J where the region where the cross section 
drops and the uncertainties grow large moves to larger 
values of prHjj or it — Acpn-jj going from CMS loose 
to ATLAS to CMS tight. We wih come back to this in 
Sec. UVBl 



2. Combination of Exclusive Cuts 

As the cases of A(j)H-jj and prHjj already show, one 
has to be careful when cutting on variables which effec- 
tively force the kinematics in the exclusive 2-jet region 
and induce large logarithms in the perturbative series. 
Whether implementing a cut-based approach or in mul- 
tivariate analysis, it is important to take into account the 
uncertainties induced by the exclusive restriction. As an 
illustration of the application of the ST method to a more 
general case, we now consider case where we combine cuts 
on both PTHjj and Acjjn-jj ■ 

Specifically, we study the exclusive 2-jet cross section 
as a function of p'rHjj with an additional constraint that 
we select only events which already have A(j)H-jj > 2.6. 
Following Eq. ([T]), the corresponding exclusive 2-jet cross 
section can be expressed as 

a2 {A^H-,j > 2.6, prHjj < pThjj ) (25) 
= cr>2 - a>3{A(f>H-jj < 2.6 or prHn > Pi^Hjj) ■ 

Taking p ~ Q for simplicity, the corresponding exclusive 
uncertainty is now given in terms of the uncertainties 
obtained by scale variation in the inclusive cross sections 
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FIG. 7: Exclusive pp ^ H + 2 jet cross section via ggF at NLO for as function of p'thh (left panels) and tt — l^cfy^ljj (right 



panels) for both ATLAS and CMS VBF selections. 



as In Fig. [SI we show a2 as a function of the p'rHjj with 

fixed A0jj_jj > 2.6 for the ATLAS and CMS loose VBF 
Al{A(t)H-j J > 2.6, pTHn < P^THjj) (26) selections. As before, the cross section for fi = uih is 

= A^^l + A^^(A0H-,,- < 2.6 or pth^ > Pt%j) ■ 
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TABLE IL Perturbative uncertainties at NLO in the exclusive pp H + 2 jet cross section via gluon fusion for cuts on pthh 
and A(f)H-jj for both ATLAS and CMS VBF selections. 




PTHn [GeV] p™^. . [GeV] 

FIG. 8: Exclusive 2-jet cross section as a function of p'rHjj with an additional cut A(f>H-]] > 2.6 using the ATLAS (left panel) 
and CMS loose (right panel) VBF selections. 



the central solid blue curve and the green dashed and 
dotted curves show the result of direct scale variation by 
a factor of two, while the outer solid orange lines show 
the uncertainties obtained from Eq. (|26p . As shown in 
Table mi for p^V^^-^. = 30 GeV we now get 56%, 31%, and 
53% uncertainty for ATLAS, CMS loose, and CMS tight, 
which is slightly increased compared to not having the 
additional cut on A(l>H-jj ■ For large values of p'rHjj 
uncertainties in Fig. [5] correctly reproduce the exclusive 
uncertainties for A2(A(/) g_jj > 2.6) without the cut on 



PTHn [see Figs. [7(b)] and [7(d)] 



B. Uncertainties in ggF-VBF Separation 



The VBF production process is characterized by two 
forward jets with large rapidity separation and large di- 
jet invariant mass. The VBF selection cuts used by the 
ATLAS and CMS experiments enhance the VBF con- 
tribution, but a significant ~ 25% ggF contribution re- 
mains. Since the VBF cross section is known rather pre- 
cisely, an important source of theoretical uncertainty in 
the extraction of the VBF signal is the large perturba- 
tive uncertainty in the ggF contribution. After subtract- 
ing the non-Higgs backgrounds (which are of course an- 
other source of uncertainty), the measured cross section 
for Higgs production after implementing the VBF selec- 
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FIG. 9: Theoretical uncertainties of the ggF contribution relative to the VBF cross section as function of p'rijjj (left panels) 
and /S.(jf^ljj (right panels) for the ATLAS VBF selection (top panels) and CMS loose VBF selection (bottom panels). The 
solid orange lines show the perturbative uncertainties in CTp'' , the green dotted lines a flat 20% parametric uncertainty in crf^^ , 
and the dashed blue lines both contributions added in quadrature. 



tion is given by 

(27) 

For the purpose of extracting the VBF cross sec- 
tion we effectively have to subtract the theory pre- 
diction for crf^(A0™_i^.^.) from cr^°^™'''='i(A(/)5J'^j.). 
Therefore, the relevant figure of merit is 
Aaf^(A<^™4.)/f^r^(A05}'4), i.e., the theory un- 
certainty in crf^^ measured relative to the expected VBF 
cross section, 0-2'^^. 

In Fig. |9] we show the ggF uncertainty relative to 
the VBF signal cross section over a range of p'thh ^'^'^ 
^(j^H-jj using the ATLAS and CMS loose VBF selections. 
In these plots, the solid orange curve shows our results 
for the NLO perturbative uncertainties (corresponding 
to the orange lines in Fig. [7]) . For comparison, the green 



dotted curve shows a fixed 20% uncertainty in the ggF 
cross section, i.e., taking Acrf^^ = 0.2 crf*^, which for 
example could be due to PDF and Us parametric uncer- 
tainties. Hence, the green dotted lines effectively track 
the size of the ggF cross section relative to the VBF cross 
section (multiplied by 0.2). In the dashed blue lines, both 
uncertainty contributions are added in quadrature. 

In the region of low p'xHjj or tt — Ac/i^J:^-^- , the relative 
uncertainty coming from the ggF contribution quickly 
increases below pxHjj ^ 30 GeV and tt — A(j)H-jj ^ 0.4. 
This is despite the fact that the relative ggF cross sec- 
tion quickly decreases there, as can be inferred from the 
decrease in the dotted green lines. In this region, the to- 
tal uncertainty shown by the blue dashed curve becomes 
completely dominated by the perturbative ggF uncer- 
tainty. Hence, one should be careful when implementing 
and optimizing either indirect restrictions on additional 
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FIG. 10: Exclusive 2-jet cross section over a range of rrij"* for fixed PTHjj < 30 GeV (left panel) and fixed A(j)H~jj > 2.6 (right 
panel) for the ATLAS VBF selection. 




FIG. 11: Perturbative uncertainties of the ggF contribution relative to the VBF cross section over a range of m^"* for fixed 
PTHjj < 30 GeV (left panel) and fixed A(l>H^j > 2.6 (right panel) for the ATLAS VBF selection. 



radiation, like A^^-jj, or explicit pT-vetoes like pTHjj, 
since the gain in sensitivity in the Higgs signal from re- 
duced non-Higgs backgrounds must be weighed against 
the increased theoretical uncertainty in separating the 
ggF and VBF contributions. 

We already saw in Sec. IIV Al that the perturbative un- 
certainties in the exclusive 2-jet cross section also depend 
on the chosen VBF cuts and increase with a higher cut on 
the dijet invariant mass, rrijj. The reason for this effect is 
that at higher rrijj the effective hard scale in the process 
is also pushed higher causing the logarithmic corrections 
at a given value of p'rHjj increase. This is seen ex- 
plicitly in Fig. [TOl which shows the exclusive 2-jet cross 
section over a range of mj^^ using the ATLAS VBF se- 
lection for a fixed cut PTHjj < 30 GeV or A(pH-jj > 2.6, 
where the curves have the same meaning as in Figs. [7] 
and m As expected, with a cut on prHjj < 30 GeV, we 
see that the relative uncertainty in the ggF cross section 
grows for larger nijj values, and reaches almost 100% for 
rrijj ^ 800 GeV. Note however that for such large rrijj 



cuts one might have to reevaluate whether /i — mn is 
still an appropriate scale choice for this process. With 
a cut on A(j>H-jj > 2.6, the relative uncertainty in the 
ggF cross section stays roughly constant for larger rrijj 
presumably because this cut is somewhat milder, which 
we also saw in the results in Table HIl 



In Fig.[TT]we show the ggF uncertainty relative to the 
VBF cross section analogous to Fig. [HI We can clearly 
see that in this case tightening the cut on rrijj does im- 
prove the separation of the ggF and VBF contributions, 
as the perturbative ggF uncertainty relative to the VBF 
cross section, shown by the orange curves, decreases. In 
this case, the overall reduction of the ggF contamina- 
tion relative to the VBF cross section is stronger than 
the increase in the perturbative uncertainties of the ggF 
contribution. 
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V. CONCLUSIONS 

In order to enhance the VBF signal over non-Higgs 
backgrounds as well as the ggF contribution, the typi- 
cal VBF selection cuts used by the ATLAS and CMS 
experiments include either indirect or direct restrictions 
on additional emissions. Such restrictions constitute a 
nontrivial jet binning, where the inclusive 2-jet cross sec- 
tion is effectively divided into an exclusive 2-jet bin and 
a remaining inclusive 3-jet bin. 

With such a jet binning one has to account for two 
sources of perturbative uncertainties. In addition to the 
absolute yield uncertainty which is correlated between 
the jet bins, there is also a migration uncertainty which is 
anticorrelated and drops out in the sum of the bins. This 
migration uncertainty is associated with the perturbative 
uncertainty in the logarithmic series that is introduced by 
the exclusive binning cut. As the binning cut becomes 
tighter, the logarithms grow large and eventually lead to 
a breakdown of fixed-order perturbation theory, at which 
point a logarithmic resummation becomes necessary. 

In practice, the experimentally relevant region typi- 
cally lies inside the transition region between the fully 
inclusive region (no binning) and the extreme exclusive 



region (very tight binning). In this region, fixed-order 
perturbation theory can still be applied. However, since 
the logarithms are already sizeable, one has to explicitly 
take into account the migration uncertainty. This can be 
achieved using the ST method. 

We studied in detail the application of the ST method 
for pp ^ H + 2 jets via ggF, including its generalization 
and validation against alternative prescriptions. We find 
that the perturbative uncertainties are very sensitive to 
the exclusive cut and can quickly become sizeable. While 
applying a strong restriction on additional emissions is 
expected to increase the sensitivity to the VBF signal, 
it is not necessarily beneficial for distinguishing the VBF 
and ggF production modes because of the quickly increas- 
ing ggF uncertainties. Hence, it would be important to 
include the perturbative uncertainties as a function of the 
binning cut when optimizing the experimental selections. 
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